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Abstract 

We study particular solutions of the inner equation associated to the 
splitting of separatrices on generalized standard maps. An exponentially 
small complete expression for their difference is obtained. We also provide 
numerical evidence that the inner equation provides quantitative informa- 
tion of the splitting of separatrices even in the case when the limit flow 
does not. 



1 Introduction 

The phenomenon of the splitting of separatrices occurs when a dynamical system 
having an invariant object (a fixed point, a periodic orbit, a torus, etc.) with 
coincident branches of its stable and unstable invariant manifolds (a separatrix), 
is perturbed. Generically, a new invariant object of the perturbed system arises 
which still possesses stable and unstable invariant manifolds but which no longer 
coincide. 

The problem of measuring the size of this splitting is long standing in Dy- 
namics. It is related to the existence of transversal homoclinic points and, 
consequently, with the non-integrability and with the size of the stochastic zone 
of the system under study. 

The most popular tool for measuring the splitting of separatrices is the Mel- 
nikov approach [Mcl63] . It is based on classical perturbation theory and provides 
a first order approximation for the splitting by using the distance between the 
stable and unstable invariant manifolds of the perturbed system. Nevertheless 
there are plenty of interesting (and in some sense generic) situations where this 
approach fails: when the Melnikov function does not predict correctly the size 
of the splitting or when no Melnikov function is available, for instance, when 
intcgrablc systems near simple resonances are perturbed. In this case, Poincare 
already detected in [Poi99] that the separatrix splits but it turns out that the 
size of this splitting is exponentially small in the perturbation parameter, what 
it is usually known as a beyond-all-orders phenomenon. Consequently a direct 
application of a first order perturbation theory never will be able to provide a 
good estimation for this exponentially small splitting. There are other settings, 
related for instance to Arnold diffusion and fluid transport, when the splitting of 
separatrices is exponentially small in the perturbation parameter, but from now 
on we will restrict ourselves to the case of near identity, analytic, area-preserving 
maps. 
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1.1 Exponentially small splitting of separatrices in ana- 
lytic maps 

Throughout this introduction wc will avoid precise statements and technicalities 
but we will give the main ideas about the exponentially small phenomenon. 

Consider an area preserving analytic map, close to the identity, that is, a 
map which can be written as 

G(z,h) = z + hg(z,h), zel 2 (1) 

where h is a small parameter and g(0, h) = 0, so that the origin is a fixed point 
for any value of h. Assume also that the origin is a weakly hyperbolic fixed 
point. Namely, redefining the parameter h if necessary, the eigenvalues A, A" 1 
of -DG(O) are of the form A = e h — l+0(h). In this case, there exist W s and W u , 
the stable and unstable invariant manifolds of the origin, respectively. The goal 
is to measure the discrepancy between these invariant manifolds. Notice that, 
since for h = the map G(z,0) = z, this is a beyond-all-orders phenomenon. 
The strategy is to not consider the first approximation of the map G as simply 
taking h = 0, but as the time h map of the vector field 

z' = g(z,0). (2) 

It can be seen, for instance [FS90] , that this approximation holds under generic 
and checkable assumptions. If the vector field (2) possesses a homoclinic con- 
nection 70 associated to the origin (the fixed point), then one expects that the 
exponentially small splitting of separatrices phenomenon arises for maps of the 
form (1). In fact in [FS90] is proved that, for any p £ W s 

dist(p, W u ) < K a e- 2 ™/' 1 (3) 

being a > and K a a constant depending on a and p, but independent of h. 
Nevertheless this upper bound is not useful for deciding whether the separatrix 
7o splits or not. It turns out to be mandatory to obtain an expression for the 
asymptotic behavior of the splitting. 

We emphasize here that, even when the distance between W s and W u seems 
a good choice for measuring the splitting, it depends on the point p. This is 
because this measure does not exploit the area preserving character of our map. 
There are several quantities more appropriate for this task. One of them is the 
Lazutkin invariant (see formula (11) in Section 2.1) which is related to the angle 
between both W s and W u at a homoclinic point. An upper bound similar to (3) 
for the Lazutkin invariant can be obtained but with K a depending only on a. 

If the asymptotic behavior for the splitting has to be proved, the first ques- 
tion that arises from (3), is how much bigger a could be. To find this optimal 
value of a one has to know the analyticity domain of 70, the homoclinic connec- 
tion of the vector field (2). It is proven in [Fon95] that 70 has complex singular- 
ities, henceforth it is analytic in a maximal complex strip {t 6 C : \hnt\ < oo}- 
The bound (3) holds for any a < do, changing K a appropriately. Notice that 
if we take any fixed er* < ero, we do not obtain a sharp upper bound, simply 
because the result also holds for a, with cr* < a < ctq and henceforth, taking h 
small enough we get a better estimate than the previous one. In consequence 
any asymptotic formula will require taking a arbitrarily close to 00 as a function 
oih. 
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The key point for proving the bound (3) is to obtain good parameterizations 
for the invariant manifolds W s , W u , which are analytic in the complex strip 
{t G C : |Imt| < cr} with a < o~$. The natural parameterization for the invariant 
manifolds are functions j n,s (t) that satisfy 

G(^ s (t),h)=Y' s (t + h). (4) 

Notice that the homoclinic connection 70 satisfies this invariancc equation for 
the time h flow of the vector field (2). As we have mentioned in the above 
paragraph, to obtain an asymptotic formula for the splitting it is necessary to 
find solutions of the invariance equation (4) defined for values of t arbitrarily 
close to do as a function of h. Since the strip is limited by the singularities of 
the homoclinic connection 70, this study becomes harder when the values of t 
are closer to these singularities. The inner equation is a suitable approximation 
of the invariance equation (4) for values of t close to these singularities. 

The main goal of this paper is to derive the inner equation for a large set of 
area-preserving maps (the so called generalized standard maps) and to obtain 
information about some special solutions and their difference. This is a first 
step in the proof of an asymptotic formula for the splitting of the invariant 
manifolds for these maps, but obtaining this formula is beyond the scope of this 
work. Nevertheless we will provide some numerical results which, combined 
with heuristic arguments (see Section 4, in particular formula (46)) support the 
relation between the splitting and the inner equation. 

1.2 The inner equation. An overview 

The study of the inner equation has been at the heart of the proof of the 
exponentially small splitting of separatrices in many examples, for maps [Gcl99, 
MSSlla, MSSllb] as well as for flows [GOS10]. 

In the case of area-preserving analytic maps, the use of the inner equation 
dates back to [Laz03], where a scheme to obtain an asymptotic formula for 
the splitting of separatrices of the Chirikov standard map was established. In 
that paper, a particular instance of the inner equation was introduced: the so 
called semi-standard map. Further development of the ideas in [Laz03] lead to 
the first rigorous proof of the asymptotic formula for the Chirikov standard map 
in [Gcl99] . A brief discussion on the splitting size on the Chirikov standard map 
can be found in [GclOOb]. From the same authors, the survey on exponentially 
small phenomena [GL01] introduces, among other things, the inner equations 
associated to polynomial standard maps, and lists in an informal way asymptotic 
formulas for the splitting of separatrices in those cases. It is also remarkable that 
in the paper [GS01] resurgence theory is applied to the study of the solutions of 
the inner equation associated to the area preserving Henon map. This paper is 
strongly related to [GclOOc]. Also in the study of perturbation of the McMillan 
map [MSSffa, MSSllb] resurgence methods were applied to study the inner 
equation. Summarizing, one can find rigorous results on the inner equation 
in [Gel99, GS01, MSSllb] in particular examples which are covered under our 
present work, which also includes and generalizes the ones present in [GL01] 
and the numerical study [GS08]. 

In the case of flows, the inner equation has been a successful tool to measure 
the splitting of separatrices when the Melnikov function fails to predict the 
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size of the splitting, like in the rapidly forced pendulum (see [GclOOa, GOS10] 
or [BFGS11] for a generalization to arbitrary polynomial Hamiltonian systems 
of one and a half degrees of freedom, following the study on the inner equation 
in [Bal06]). A different technique based on continuous averaging to study the 
exponentially small behavior of the splitting can be found in [Trc97]. 

The purpose of the present paper is twofold, a combination of rigorous the- 
oretical results in a general setting and numerical experiments avoiding lengthy 
proofs in particular examples. One of the numerical examples shows a a type 
of behavior that is not covered by the surveys [GL01, GS08] (see the end of this 
section). 

We study some second order difference equations, called inner equations, 
which have the form either 

4>{z + 1) - 2<j>{z) + <j>(z - 1) = -<j> n (z) + G{^z)) 

or 

<f>{z + 1) - 2<jy(z) + <j){z - 1) = -e"* (z) + G(e*( z )), 

depending on the class of maps under consideration, and where G(w) is an 
analytic function such that G(w) = 0(w n+1 ). 

These equations appear, in particular, in the problem of exponentially small 
splitting of separatriccs in generalized standard maps (see next Section for defi- 
nitions) , but they can appear in studies of other types of maps (with parabolic 
fixed points, for instance), and, with this applicability in mind, we consider 
them in their full generality (see equations (27) and (28)). In particular, our 
present results generalize those on the inner equations appearing in [Laz03, 
Gel99, GelOOb, GL01, GS01, MSSlla, MSSllb]. It is important to remark 
that in the previous literature on the subject the symmetries of the particular 
problems under consideration where exploited extensively in the proofs. Our 
present formulation does not rely on additional symmetries, making it suitable 
for applications. In particular, wc provide all the technical details and complete 
proofs of the statements concerning the inner equations and their solutions. As 
a side comment for the specialists, there are several technical improvements in 
the proofs of our theoretical results, which we expect can be applied in other 
problems related to difference equations. 

We describe a large set of formal solutions of these inner equations, from 
which some true solutions are obtained, and we derive a complete formula for 
their difference. The main results are collected in Section 3, while Section 2 plays 
the role of a more detailed introduction of the problem and description of some 
of the known results. Sections 5, 6 and 7 are devoted to prove the theoretical re- 
sults, while Section 4 contains the numerical results with a non-rigorous exposi- 
tion of their relation to the developed theory. It should be remarked that the re- 
lation between the inner equation and the actual computation of the splitting, in 
the particular cases where proofs are available (see [Gel99, MSSlla, MSSllb]), 
is lengthy and full of technicalities. Our exposition here tries to give the reader 
an idea of the link between the inner equation and splitting size, by making very 
strong assumptions, in order to explain the obtained numerical results. These 
assumptions are fully proved in the literature for the Chirikov standard map 
and the McMillan map. 

The numerical experiments have been conducted to test the applicability of 
the theoretical results. Although academic in nature, they show the relation 
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between the splitting of separatrices and the difference between two solutions 
of the inner equation. Moreover, the main example exhibits a behavior that is 
not covered by the surveys [GL01, GS08]. In this example, given by the map 

x\ fx + y + e(x — x 3 ) — e 2 x"' 
yj I y + e(x — x 3 ) — e 2 x 7 

where e is a small parameter, although the size of the splitting is much larger 
than the guess suggested by [FS90] , the leading term of its asymptotic behavior 
is provided by the inner equation. As a matter of fact, the splitting size in this 
example behaves asymptotically when s — > as 

A ( 7T 2 2 5 /40fr(3/4) 2 )\ „ , . , 

-exp — — H —j- (1 + higher order terms), 



h™l 3 ' \ h ' h 1 / 2 

where e = 4 sinh 2 (/i/2) and A is a constant related to some inner equation, while 
the naive guess provided by the limit flow (see Section 2 for details) , in this case 
the Duffing equation x = x — x 3 , would be exponential with exponent — n 2 /h. 
That is, the correction term is larger than any power of h. See Sections 2.4 
and 4. 

We remark that although the computation of the actual splitting has been 
performed by using the multiple precision package PARI-GP, the computation of 
the leading term has been achieved by using the standard long double precision 
in C. 



2 Generalized standard maps and exponentially 
small splitting of separatrices 

2.1 Generalized standard maps 

We will say that an area preserving map (x*,y*) = F{x,y) is a generalized 
standard map if it can be written in the form 

x* = x + y + f(x,h), ^ 

y* = y + f(x,h), 

where ft, is a small parameter. We will assume that / depends analytically in its 
arguments in \h\ < ho, \x\ < po, for some fixed ft-Oi Po > 0. We will be interested 
in the case when the origin is a fixed point of F, that is, /(0, h) = 0. Moreover, 
we will assume the origin to be weakly hyperbolic, although our study may be 
applied also to the case of a parabolic fixed point. 

The parameter h is chosen in such a way that spec£>F(0,0) = {e h ,c~ h }. 
This last condition is equivalent to impose f'(0,h) = ^/(0, h) = e, with e = 
4sinh 2 (/i/2). We further assume that 

f(x, h)=Y / fk(x)h k+2 = efo(x) + 0(h 3 x). (6) 

Under these conditions, the map (5) can be written as a close to the identity 
map: with the scaling x = x, hy = y, it becomes (using again x and y as 
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variables) 

j x* = x + hy + 0(h 2 x), 

\y* = y + hf (x) + O(h 2 x). 

When h is small, the map (7) is well approximated by the time h map of the 
flow of the Hamiltonian system 

I x = II. 

(8) 




We assume that the origin in (8), which is a fixed point, possesses a homoclinic 
connection, 7o(t) = (xo(t),yo(t)). By a shift on t, we can choose 70 such that 
xq is an even function, that is, 70 intersects transversally the line {y = 0} at 
t = 0. The invariant manifolds of the origin for the map (7) are close to this 
homoclinic connection. Hence, if h is small, by the conservation of the area, 
they must intersect. It is not difficult to check that the expansions in powers of 
h of the stable and unstable curves coincide. As a consequence, the expansion of 
the angle of intersection in powers of h vanishes, which, in view of the analytic 
nature of the problem, suggests that this angle may have an exponentially small 
behavior in h. In fact, Fontich and Simo, in [FS90], obtained an exponentially 
small upper bound for the angle. They showed that if 70 is analytic in the 
complex strip {|Imi| < <7q} and the map F is defined around the homoclinic 
orbit, then, for any < a < ctq, the distance between the stable and the unstable 
manifold of the origin of (7) is bounded by K a e- 2 ™/ h , for any < h < h a , 
where K a and h a are positive constants depending on a and K a depends also on 
the point where this distance is measured. Restoring to the original variables, 
the same applies to the invariant manifolds of the origin of (5). 

Equivalently, a natural parametrization j(t) = (x(t),y(t)) of the invari- 
ant manifolds of the origin of (5), when condition (6) is satisfied, that is, a 
parametrization satisfying F o j(t) — -f(t + h), must be a solution of the differ- 
ence equation 

x(t + h)- 2x(t) + x(t -h) = f(x(t),h), (9) 

with y(t) = x(t) — x(t — h). This equation implies that the curve 7 = (x,y) 
is invariant by F and the action of f on 7 is conjugated to the shift on the 
parameter t: t n- 1 + h . One must supply additional conditions on 7 to obtain 
the invariant stable and unstable curves: if 7 is the unstable manifold (resp. 
stable) of the origin then limt-y-oo x(t) = is required (resp. lim^oo x(t) = 0). 
Since the left hand side of the invariance equation (9) is formally 

x(t + h) - 2x(t) + x(t - h) =4sinh 2 (~) (x)(t) = h 2 x{t) + 0(h 4 ), 



2dt, 

it can be approximated, when h is small, by the second order differential equa- 
tion 

i = fo(x), (10) 

which is nothing more than (8). 

In order to measure the difference between the invariant manifolds, it is often 
used the Lazutkin invariant at a homoclinic point p = 7"(0) = 7 s (0), 

.( P) =det(| 7 «(0) ! |7 s (0)) j (ID 
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being j u ' 3 (t) natural parametrizations of the unstable and stable manifolds. 
Unlike the angle between the invariant curves, w(p) is a symplectic invariant 
and only depends on the homoclinic orbit, not on the specific point p. Another 
symplectic invariant quantity that can be used to measure the splitting of the 
separatrices is the area of the lobe between two consecutive homoclinic points. 

Since an upper bound of the splitting of the separatrices is known, the ques- 
tion of its asymptotic behavior when h tends to arises. Some well known 
examples in the literature where this formula is available arc briefly summa- 
rized in the next subsection. 

2.2 Examples of generalized standard maps with exponen- 
tially small splitting of separatrices 

There are not many examples with a complete proof of an asymptotic formula for 
the splitting of separatrices in area preserving maps. Here we quote two. There 
is a more abundant literature about splitting of separatrices in Hamiltonian 
systems with one and a half degrees of freedom (see [SMH91, DS97, DGJS97, 



The first example is the Chirikov standard map, introduced by Chirikov as 
a basic model of the motion of a system close to a nonlinear resonance (see, 
for instance, [Chi79]). It corresponds to take f(x,h) = esin(a;), with e = 
4sinh 2 (/i/2). This map is in fact defined in the annulus, and the limit flow (8) 
is a pendulum with the saddle at the origin. The separatrix of the pendulum 
is analytic in the strip {|Imi| < tt/2} and has a singularity at t = in/2. The 
symmetries of the problem imply that there is a homoclinic point p on the line 

x = 7T. 

In [Gcl99], Gelfreich proved, following the scheme developed by Lazutkin 
in [Laz03], that 



where the series in the right hand side is asymptotic. In particular, the exponent 
in the exponential is well predicted by Fontich-Simo theorem in [FS90]. 

The second example is the perturbed McMillan map. The McMillan map 
itself was introduced in [McM71] in the modelization of particle accelerator 
dynamics. In [DRR98, MSSlla, MSSllb], perturbations of the McMillan family 
of the form 



are considered, with V(y) = J2k>2 Vky 2k analytic in a neighborhood of y = 0. 
In the above formula, h is the Lyapunov exponent of the origin, which is the 
small parameter, and e is independent of h and not necessarily small. The 
McMillan map is obtained when s = and is intcgrable with a polynomial first 
integral. See [DRR98] for more details about the McMillan map. 

With a linear change of coordinates, the map (12) can be written in the 



Tre97, LMS03, OSS03, DG04]) 





(12) 
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form (5) with 

= e(x - (2 - 4eT/ 2 )x 3 ) - e 2 (x 2 - (2 + 6eF 3 )a: 5 ) + 0(e 4 ), (13) 

where, again, e = 4sinh 2 (/i/2). The limit flow (10) is the Duffing equation 

x=x - (2-4ey 2 )x 3 , 

with homoclinic x (t) = a/ eosh(i), a = [X—leV^ 1 ! 2 (assuming |e| < (2y 2 ) _1 )- 
Its closest to the real line singularities are located at ±i7r/2. In [MSSlla, 
MSSllb], improving a partial result in [DRR98], it was proven that, if V(27r) ^ 
0, where 

is the Borel transform of V, then the invariant manifolds to the origin of (12) 
split when i ^ and the Lazutkin invariant of a particular homoclinic orbit 
satisfies 

1 k>0 

where the functions are analytic around e — 0, /3 2 = 1 — 2eV2/cosh/i 
and Bq(e) = 47r 2 V^(27r) + 0(e). If the map is written in the form (5), with 
the function / given in (13), the Lazutkin invariant has an additional h 2 in 
the denominator. Again, the exponent of the exponential is well predicted by 
Fontich-Simo theorem. 



2.3 Numerical studies for polynomial generalized stan- 
dard maps 

In [GS08], Gelfreich and Simo presented a detailed numerical study of the 
splitting of the separatrices of the generalized standard map (5) in the case 
f(x,h) = ep(x), with p{x) — 'Y^=iPkX k & polynomial of degree n with p\ = 1 
(which implies f'(0,h) = e) and p n < 0. Is is also assumed that there is a 
homoclinic curve to the origin in the limit flow system (10). 

Then, via numerical experiments, the authors showed that the asymptotic 
behavior of the Lazutkin invariant depends only on the relative position of the 
singularities of the homoclinic solution of (10), on the degree n of the polynomial 
p and on the coefficient p n : 



\Pn\ V ' 2 h» 

where v = 2(n + l)/(n — 1), p is the minimum distance to the real line of the 
singularities of the homoclinic of (10), u>(h) ^ is either a constant, a periodic 
function or a quasiperiodic function of 1/h, depending only on the number of 
singularities at |Imt| = p and their relative positions and C n depends only on 
n. 

Also in this case, the exponential behavior is well predicted by Fontich-Simo's 
theorem. 
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2.4 A discrepant example. Numerical observations 

We introduce the generalized standard map (5) induced by 

f{x,h) = s{x-x 3 )-e 2 x 7 . (14) 

Note that this map possesses terms in e 2 , like the McMillan map has (see (13)). 
Unlike the McMillan case, the function defining this map is entire. 

The limit flow (10) for this map is also a Duffing equation, in this case 
x = x — x 3 , with homoclinic xo(t) = \/2/cosh(t), whose singularities are located 
at the same place of the homoclinic of the McMillan map, being ir/2 their 
minimum distance to the real line. Hence, one could be tempted to infer that 
the exponential behavior of the Lazutkin invariant is of order e _7r l h . 

However, our numerical experiments suggest that the Lazutkin invariant at 
the first homoclinic point over the line y = 0, in the topology of the unstable 
manifold, behaves like 

e -2^)A + ... ; (15) 



/jlO/3 



where 



2 */ir 



and A = 871.683 .... In particular, the size of the Lazutkin invariant is much 
larger than the naive guess, which, in turn, suggests that the approximation of 
the invariant manifolds provided by the limit flow (10) is not good enough to 
predict the asymptotic formula of the splitting. Section 4 is devoted to explain 
these numerical experiments. In particular, we will conjecture the source of the 
function p(h) and the origin and computation of the constant A. 



2.5 Inner equation for generalized standard maps 

In all the aforementioned examples, the constants u>q, B^(e), C n and A in the 
leading term of the asymptotic behavior of the Lazutkin invariant are related 
to a suitable inner equation, whose solutions provide better approximations 
of the invariant manifolds for values of t in some regions of C than the one 
provided by the limit flow (10). Even in the case of the generalized standard map 
defined by (14), where the limit flow (10) does not provide enough information, 
the numerically evaluated constant A in (15) is obtained from such an inner 
equation. 

In order to be able to construct the inner equation we will impose several 
conditions to the function defining the generalized standard map. 

Let F be a generalized standard map of the form (5), induced by a func- 
tion f(x,h) = J2 k>a fk(x)h k+2 , satisfying the hypotheses in Section 2.1. We 
furthermore assume: 

HP1 For each k > 0, f k (x) = fk,jX j , with / Mfc + 0. 

HP2 The function k M> (d k - l)/(k + 2) has a global maximum on N. Let / C N 
be the set where this maximum is achieved. 
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Hypothesis HP2 implies a restriction in the rate of growth of the degree of 
each of the polynomials which can be at most linear in k. We also remark 
that, combining hypothesis HP1 and HP2 with the fact that / is analytic in 
the bidisk H> po x D^,,, one obtains that the domain of analyticity with respect 
to x depends on h and tends to be the whole complex plane when h tends to 0. 

We fix x G C. We introduce the new unknown <fi(z) defined by x(x + hz) = 
h- a X(j)(z), with 

k + 2 * , T 

a = — , tor any k €E 1 

dk - 1 

and A a parameter to be determined later. Note that, by definition of / in HP 2, 
a is indeed independent of k G I. The invariance equation (9) becomes 

<j){z + 1) - 2<f>(z) + <j>{z - 1) = h a \- 1 f{h- a \(t>{z), h). (17) 

With the standing hypotheses, the right hand side above admits an expansion 
of suitable positive powers of h as follows: 

fe>0j=l 

(d h -i 
3 = 1 

= h, dk \ dk ^^ dk + 0(/i rain{1 ' Q} ), 
kei 

where in the last equality we have used the definition of a and I. The inner 
equation is obtained by keeping only the first term in h in the right side of (17): 

<f>(z + 1) - 24>(z) + <)>{z - 1) = Y, /M fe ^~V dfe . (18) 

fee/ 

Let n = min{<ifc : k € I}. To simplify the notation we introduce the coefficients 
Gk such that 

kei k>n 

Now we take A such that A n_1 = — (G n ) _1 . With this choice, the inner equation 
associated to the generalized standard map is 

cP(z + l)-2<p(z) + <l>(z-l) = -cj> n (z)+ Y G k<f> k (z) (19) 

k>n+l 

being Gk = GfeA fe_1 . Notice that G{4>) := J2k>n+i Gk<j> h is analytic in a neigh- 
borhood of <f> = 0. 

In the trigonometric case one can proceed analogously. Indeed, assume that 
f(x, h) = Y.k>n fk{x)h k+2 , with / satisfying: 

HT1 For each k > 0, fk(x) = Y^j=-d k Sh,j e ^ x is a trigonometric polynomial of 
degree d k > a, with f h>dh ^ 0. 
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HT2 The function k i-> dk/(k + 2) has a global maximum on N. Let / C N be 
the set where this maximum is achieved. 

For any x € C, we define <p(z) by x(x + hz) = — ilog (h a X) + i(j){z) with 

k + 2 

a = — — , for any fc e I (20) 

and A a parameter. Then, the invariance equation (9) becomes 

(j)(z + 1) - 20(z) + (j)(z - 1) = -i/(-ilog (/i Q A) + i0(z), /i). (21) 

As in (18), the inner equation is the above equation when h — > 0. In this case, 
taking A appropriately, one obtains 

4>{z + 1) - 2cf>{z) + 4>{z - 1) = -e^W*) + G k e k ^ z) , (22) 

k>n 

where n — 1 = minjdfc : fc G /}. The discrepancy in the definition of n in both 
cases allows us to make an unified treatment of the problem in next sections. 

Since the original invariance equation (9) is autonomous, the inner equa- 
tion (19) and (22) does not depend on the choice of the complex number \ 
introduced with the new unknown cf>. Nevertheless, this complex number is es- 
sential when the size of the splitting of separatrices is studied and has to be 
well chosen. Roughly speaking, it will measure the exponentially smallness of 
the splitting which turns out to be 0{h y e- 27!lrax / h )) for some v £ R. This 
asymptotic behavior has only been proved for particular maps (see Section 2.2) 
but there are numerical evidences, Sections 2.3 and 2.4, that it also holds in a 
more general setting. We plan, in a future work, to prove it for the generalized 
standard maps. 

In the examples presented in Sections 2.2 and 2.3, \ is chosen to be the 
location of the closest to the real line singularity of the homoclinic solution 
70 of the limit flow (10). In example in Section 2.4, x i s a l so related to the 
singularities of a homoclinic solution of some flow, which is not longer (10) but 
x = x — a; 3 — ex 1 . 

2.5.1 Some examples of the inner equation 

Here we show how the inner equation is derived for some examples. 

The first one is the map introduced in Section 2.4. Its inner equation is 

(Xz + 1) - 20(z) + 0(2 - 1) = -0 7 (z). (23) 

Indeed, in this case f(x, h) = s(x — x 3 )+e 2 x 7 with e = 4sinh 2 (/i/2) . Therefore, 
fa(x) = x - x 3 , f2k(x) = f2k,ix - f 2 k,3X 3 - !2kW and f2k-i(x) = 0, for k > 1, 
which implies that do = 3, dik = 7 and cfefc-i = for k > 1. In this situation, it 
is clear that n = 7, a = 2/3 and the set I = {2}, therefore the right hand side of 
equation (18) is /2.7A 6 7 and defining A adequately we encounter equation (23). 

Now we compute the inner equation for the generalized standard map in- 
duced by f(x, h) = e(x - x 3 ). In this case fak(x) = f2k,\x - f2k,3% 3 , d 2 k = 3, 
f2k+i(x) = and c^fc+i = 0, for k > and this implies that n = 3, a = 1 and 
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the set I = {0}. Then, the right hand side of equation (18) is /o,d ^ 2( / )3 an d we 
obtain the inner equation 

4>{z + 1) - 24>{z) + 4>{z - 1) = -(f> 3 (z). (24) 

We can also encounter inner equations having infinite terms in its right hand 
side. For instance by considering f(x, h) = esin(x) + £fc>i akh 2k+2 sin ((fc + 
l)x). In this case d^k = 2fc + 2, d-2k+i = 0, n = 2, a = 2 and I = {k E N : 
k is even}, so that the inner equation is 

(j>{z + 1) - 24>{z) + <\>{z - 1) = -c^ + G k (c k ^). 

k>2 

The main purpose of this paper is to provide some particular solutions of 
the inner equation (19) and (22) as well as to compute an explicit formula 
for their difference. The precise statement is placed in next section, while its 
proof is spread along the subsequent ones. As we have already commented in 
Section 1.2, this computation has been in the heart of the proof of the splitting 
of separatrices in all the known examples, and it also gives an explanation to 
the numerical results concerning the example in Section 2.4. 

3 Main results 

We consider the linear operators 

&{4>){z) = 4>{z+l)-<f>{z) (25) 

and 

A 2 (</))(z) = A(0)(z) - A(0)(* - 1) = ^(z + 1) - 24>{z) + 0(z - 1) (26) 

and two types of inner equation. The first one, under the hypotheses HP1, 
HP2, which from now on we will call polynomial case, 

A 2 (</>)= ff (</>, M ) :=-4> n + G(<l>,n) (27) 

and the second one, under the hypotheses HT1, HT2, which we will call 
trigonometric case, 

A 2 (0) - 5 (0, fi) := _ c ("-^ + G(e* ft), (28) 

with G an analytic function in some open bidisk D(g) x D(/io) S C 2 and such 
that 

G(y,fi)= ^ Gk{fi)y k in the polynomial case (29) 

fe>n+l 

G(e v ,fi) = Gk(fi)e ky in the trigonometric case. (30) 

k>n 

The parameter ft is included for the sake of completeness and it is a regular 
parameter. 
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Remark 3.1. Let aeKbe such that an > 1. If we consider inner equations 
of the form either A 2 (</>) = g((f> a , /x) in the polynomial case or A 2 (4>) = g(a(j>, fj,) 
in the trigonometric one, the results in this section also hold true with the same 
proof. However, in order to avoid a new parameter, we restrict ourselves to the 
hypotheses above. 

In this section we present the results dealing with both formal and analytic 
solutions of the inner equation. 
Given v > 0, we will denote by 

c[[*-1, M] = | 00*) = E CJ ^0^ I : -+ c 
[ fc>i z 

the space of formal power series in z~ v without constant term, whose coefficients 
Cfc_i, depend analytically on \i 6 B(^iq). 

Proposition 3.2. Let n > 2, r = 2/(n - 1). 

1. If n is even, the equations (27) and (28) admit a unique formal solution 
<j> such that <j> € C[[z _r ], {fi}] with Cq 1-1 = — r(r + 1), in the case of (27), 
and, in the case of (28), <f> — 4>q g C[[z _r ], {/_*}], tuii/i 

0ow=^ T iogf — ^y (si) 

71—1 \ 77, — 1 Z^/ 

Moreover, any formal solution of the inner equation (27) belonging to 
C[[z~ r / 2 ],{ii}] is of the form (j>(z — c, fx) , for some c £ C and cq such that 
Cq -1 = — r(r + l). TTie same applies to any formal solution </> o/ (28) sitc/7 
that<j>-4> € C[[z~ r /%{(!}]. 

2. If n = 2m — 1 wif/i m > 2 , £/ie formal solutions are 

^'^ = E^ a ,.,://; In-' c. (32) 

^ 0<j<[^] 

in i/ie case of (27), wii/i Co = Cq.q satisfying Cq -1 = — r(r + 1), and 

^( z ^) = ^ T iog(-^- T ^) J2 cfc_ij(/i)W«, 

fc ^ 0<j<[^] 

(33) 

in i/ie case o/ (28). TTie symbol [x] stands for the integer part of x. The 
coefficients Ck-i,j are analytic functions in B([1q). 

The solution is unique provided that c m _i.o = 0. Any other formal solution 
of the form (32) or (33) is obtained from these ones by translation. 

Now we deal with the analytic solutions of the inner equation. Let us define 
the complex domains where these solutions are defined. For any p, 7 > 0, we 
introduce 

D* >p = {zeC: |Imz| > - 7 Rcz + p}, D» p = -D s J>p . (34) 



Inner equation for area preserving maps 



14 



D u 
i-p 


Imz 

Rez 


7 -ip 





Figure 1: Unstable domain 



Let (f>o be defined by (31) in the trigonometric case and <po = in the 
polynomial case. Let 4>q be the truncation up to order n in z~ r of the formal 
solution provided by Proposition 3.2, that is, if n = 2m with m > 1, 

n 

0o(2) = 0o(^) + X] C fc-l(M>" fcr ' ( 35 ) 
fc=l 

and, if n = 2m — 1 with m > 2, <f>o in the polynomial case is 

" 1 

<M*0 = J! 37 H Ck-iAf*) lo S J z ( 36 ) 
fe=lZ o<i<[i^] 

and in the trigonometric case 

0o(z) = 0o(*)+X,^ ]T c fe _ 1 . J ( M )log''z. (37) 

Theorem 3.3. (Existence theorem) Let r = 2/(n — 1) and Co fee smc/i i/iai 
Cq _1 = — r(r + 1). -For an?/ 7 > there exists po big enough such that for 
any p > po, the inner equations (27) and (28) have two analytic solutions 

u ' s (z,/i) = ^o(2) + ^ u ' s (^^), 

SUp |z r+2 -0 U ' 8 (z,/i)| < +OO, 

Now we state the theorem for the difference <fi n — <fi s . First we define the 
complex domain (see Figure 2) 

E^ p = D^ p r\D s l p n{z g C : Imz < 0}\{z e C : |Rcz| < 1, |Imz| < p+j}, (38) 

where the difference between two solutions of the inner equation (27) 4> u — cf> s is 
defined. 
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Im 2 



Figure 2: Inner Domain 



To unify the notation we introduce the new parameters 



1 = 



r + 2 polynomial case 
2 trigonometric case 



, Co polynomial case 
1 trigonometric case 



Theorem 3.4. Let cf> n ' s be two analytic solutions of equations (27) and (28) 
satisfying the conditions stated in Theorem 3.3. 

Their difference <p u — <jf : E^ jP x -B(^o) C, can be expressed as 

r(z^) ~ ^{z^) = C 1 (z,^P 1 M^ lkz + C2(z,^P 2 M^ lkz (40) 

fe<0 fe<0 

with p\,p\ analytic functions in B(^lq) and (11C2 satisfying that: 
1. their wronskian 



(l(z,p.) (2(z,[l) 

d(z + l,p) ( 2 (z + l 1 ri) 



= 1 



2. there exists a constant C such that for any z € Ej_ p and jj, £ B(p,o), 



z^e^id&u.yd^iz,^ <C, 



los 



— Ca(2,/i) 



rd t {2l- 1) 



< C 



wii/i v = I — r , er = if n > 3, u = t — 1 if n < 3 ; (7 = if n = 2 and 
(7 = 1 if n = 3. 

From now on we will skip the dependence on being always analytic. 



4 Numerical results 



In this section we present some numerical results concerning the generalized 
standard map (5) given by the functions fi(x, h) = e(x — x 3 ) — e 2 x 7 in (14) and 
/ 2 (x, h) = e(x - x 3 ). We recall here that e = 4 sinh 2 (/i/2). 
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We notice that both functions /i,/2 satisfy the hypotheses of Section 2.5. 
Henceforth, as we show in Section 2.5.1, we can construct the inner equation 
for the generalized standard map induced by them: 

A 2 (0) = ~<j) 7 and A 2 (0) = -0 3 . (41) 

The first one corresponds to f\ and the second one to f 2 - 
Let 

e - 4> u - 4> s , 

be the difference between the two solutions of the inner equation (41) given by 
Theorem 3.3. First of all, in a general setting, we relate the main term of 
with the Lazutkin invariant for the standard map (5) induced by /. Next, we 
compute the actual Lazutkin invariant for the maps defined by /i and f 2 which is 
computed numerically by using multiprecision routines. After we summarize the 
method to compute the main term of the difference O := U — <jf , by exploiting 
the theoretical framework we have developed. One aspect worth to remark is 
that these computations have been performed through standard long double 
precision arithmetic. 

A similar, but more detailed, numerical comparison between the Lazutkin 
invariant and the difference is performed in [GG11] for the Swift-Hohenbcrg 
equation. 



4.1 The relation between the Lazutkin invariant and 

For computing the first asymptotic term of we now take advantage from the 
fact that we have an alternative expression for by using the functions £i and 
(2 given in Theorem 3.4. Indeed, we actually can write the difference as 



0(z) = Ci(z) ^Pie 2 * + C2(z)J2pI c2 



k<0 k<0 



with Pj(z) = J2k<oPk e2 ^ lkz ' 3 = 1)2 1-periodic functions. We recall that by 
Theorem 3.4 C2) = 1 and henceforth Pl = W(Q, C2) and p 2 = W(&,Q). 

On the one hand, we introduce the new quantity uji n (z): 

Uin{z) := ~Y z W{Q ^ l){z) = T z P2{z) = E 27ri ^ e2 "" Z « -2^ 2 -ie- 27riz . 

k<0 

(42) 

The last equality has been deduced as Imz — > —00. On the other hand, note 
that by using the first approximations of (1 and (2 in Theorem 3.4, since £ > 0, 
Ci(^) — > as Imz — > —00 and (2(2) = 0(z ), the main term of is 

0(Z) = C1(Z)P1(Z) + Q 2 {Z) P 2{Z) « C2{Z)P2{Z) « Z \ d ^_ x f 2 ™- 

We recall here that only p 2 _ 1 is unknown, the other quantities are defined in 
terms of the inner equation. Henceforth, both uJi n (z) and z~ £ Q(z), are asymp- 
totically equivalent. 

In order to compare the numerical results with our theoretical framework 
we will gather in a rather informal way several facts, some of them not proven. 
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In particular, to transform assumptions (Al) and (A2) below into proven facts 
would require involved arguments even for particular cases. By this reason, we 
will avoid precise statements. The chain of reasoning is a slight modification of 
the one in [MSSlla], which also follows [Laz03, Gcl99]. 

Let / be a real analytic function satisfying the hypotheses in Section 2.1 and 
Section 2.5. We first remark that there exists a solution of the invariance equa- 
tion (9) induced by /, x u (t), i7r-antiperiodic, entire and real analytic in t, such 
that lim Ret _>_ 00 ;E u (i) = and x u (0) = x u {-h) (and x u (t) - x u (t - h) > 0, 
for t < 0). Then, the function x s (t) = x u (—t) is also a solution of (9), 
with the same regularity, satisfying linract-yoo x s (t) = 0. Hence, 7 u,s (i) = 
(x u ' s (t), x n,s (t) — x u ' s (t — h)) are natural paramctrizations of the invariant man- 
ifolds of the origin. We notice that p = 7 U (0) = 7 S (0) = (x u ' s (0),0) is the first 
homoclinic point. Let D(i) = x s (t) — x u (t). 

Using the h-step Wronskian 



W h (u,v){t) = 



u(t) v(t) 
u{t) - u(t - h) v(t) - v(t - h) 



u(t) v(t) 
A h u(t) A h v(t) 



the Lazutkin invariant (11) can be written as 



w(p) - det( 7 ",r)| t=0 = ldet( 7 ",7 s -7")| t =o = jW h {x\D) ]t=[) . (43) 

Since both x u and x 8 are solutions of the second order difference equation (9), 
their difference D also satisfies a linear second order equation, namely 



s)x u (t),h)ds) D(t). (44) 



Notice that, if x u is close to x s , then equation (44) is close to the linearization 
of the invariance equation (9) around x n . Hence, our first assumption is that 

(Al) there is a (real analytic) solution r/i of equation (44) close to x u . 

Let r]2 be another (real analytic) solution of (44) with Whifji,^) = 1, which 
can be obtained by the "variation of constants" method. Hence, we can write 
D = c.\f\\ + C2V2 where c\ and C2 are the /i-periodic functions c\ = Wh{D, 772) 
and C2 = Wh{r}i,D). Substituting this expression for D into (43) and using that 
771 is close to x u we have that 

w(p)w^W h (f7i,D)| t=0 . (45) 

Since / satisfies the hypotheses of Section 2.5, we can construct an inner 
equation associated to the standard map induced by/. The second assumption 
is 

(A2) there exists x •= C (which can depend on h) such that for values of t 
satisfying \t - x\ = 0(h), x u ' s (t) are close to h- a X(/> u ' s ((t - x)/h). Here 
<fi n,s are the solutions of the inner equation (27) given by Theorem 3.3 
and a, A are both parameters introduced in Section 2.5. Since / is real 
analytic one can assume that Im% > 0. 
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As a consequence, since, by Theorem 3.4, (i(z) = d z (f> u (z) + 0(z r+1 e 27Tlz ), 
± u (i) « h-^X^dt - X )/h) « h-^Xddt - X )/h). 

Recall now that pi(z) = —W(Q,Ci)(z). Hence, taking into account the scaling 
and assumption (Al), for values of t close to x, 

W h ( m ,D)(t) « W h (h- a - 1 \( 1 ,h- a \0)((t - x )/h) = h- 2a - 1 X 2 P2 {{t-x)lh). 

Then, since Wh{r]i, D)(t) and W((i, Q)((t — x)/h) are both /i-periodic and that 
the first one is a real analytic function we easily have that, for real t 

j t W h {x«, D)(t) a 2h~ 2a ' 2 Rc (V • ±p 2 ((t - = 2h~ 2a - 2 Re (A 2 • Win ((i - X )/h)) 

with Wj n defined in (42). Hence, evaluating at t = 0, 

w(p) w 2h- 2a - 2 Re (A 2 • Win(-x//i)) ■ (46) 

Our goal now is to check numerically the above formula for the maps induced 
by h an d h- 



4.2 The limit flow and its singularities 

In the cases of the Chirikov standard map and the perturbations of the McMillan 
map in [Gcl99] and [MSSlla], resp, x = iT/2 is the closest to the real line 
singularity of the homoclinic orbit of the limit flow (10). In the maps induced 
by fi(x,h) = e(x — x 3 ) — e 2 x 7 and f2(x,h) = e(x — x 3 ) under consideration, the 
closest to the real line singularity of the homoclinic of the limit flow x = x — x 3 
is also iTr/2 (see Section 2.4). Nevertheless, our numerical computations show 
that it is not the right guess for x in the case of /i . For this reason, we consider 
the higher order (in h) limit flow 

x = x - x 3 - ex 7 . (47) 

The parametrization, xo(t,h), of the homoclinic loop to the origin such that 
.To(0, h) = has a singularity at 

dx 

p(h) 



x o (o,/0 y/x 2 /2~xi/4-sx s /8 

where xo(0, h) = \[2 + 0(h 2 ) is the positive root of x 2 /2 — x 4 /4 — ex 8 /8 and the 
integral is computed along the real line. The other singularities can be obtained 
changing the path of integration. It can be seen that 

p{h) = i| - i 2V4l W 4 ) 2 ftV2 + {hV 2 ). (48) 
2 v 71 " 

We remark that, although the singularities of the homoclinic of (47) tend to the 
singularities of the limit flow 'x = x — x 3 (10) (in a rather slow way), they are of 
a different type: whereas the latter are poles, the former are branching points. 

We choose the values x — if — i- — ^j/ 4 ^ 1 fe 1 / 2 for f 1 and x = if f° r fi an d 
we will assume that (i) holds for them. 
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4.3 Numerical computations 

Wc define now 

w(h) = h 2a+2 \- 2 c 2 ^/ h Lu(p), Qiniz) = 2e 2Ti2 Re (w in (z)) (49) 

taking A = 1 and on the one hand a = 2/3 for fi and on the other hand 
a = 1 for f-2- We note that, since \ nas n0 rea l part, checking formula (46) is 
equivalent to check that 

ui(h) ~ uii n ( — x/h) lim u>(h) — lim uj- ln (z). 

First we show the results for co(h). We have computed numerically this 
quantity by using multiprecision routines written in PARI-GP. In the following 
figure we show the computed values for f±(x,h) = e(x — x 3 ) — e 2 x 7 and for 
the map induced by f2(x, h) = e(x — x 3 ). Let us denote by Qi(h) the value of 
Q(h) for the corresponding maps fi, i — 1,2. We have added a correction factor 
esc = 85 ■ 10~ 4 in order to have the same magnitude for both values of Q(h). 



6.755 




6.725^ , , , , ; 

0.0005 0.0015 0.0025 0.0035 0.0045 0.0055 

h 



These numbers have been obtained computing explicitly uj(p) = dct(7"(0), 7 S (0)), 
following the strategy in [DR.R99]. Due to the exponentially small behavior of 
this quantity, it has been necessary to compute 7"' s (0) with increasing accuracy, 
thus making impossible to achieve very small values of h. 

Notice that, in the case of the map induced by f2(x, h) = s(x—x 3 ), the values 
of 002(h) converge quite fast, when h becomes smaller, to a constant value 

uj 2 {h) w 1.00083- 10 5 (50) 

In the case of the map induced by f\(x, h) = e(x — x 3 ) — £ 2 x , the convergence of 
the values oiw\(h) is slower, as the previous figure shows. However, computing 
uj\(h) for h = l/2000 + A:/40000, k = . . . 199, and making some assumptions on 
the form of the asymptotic expansion of u>i(h) in h, it is possible to extrapolate 
the limit value with better accuracy. 
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In this way, we have obtained that it is 

wiO) w 871.683 



(51) 



We remark that, with the computed data, in which each value of Qi(h), i = 
1,2, has a few hundreds of correct digits, it would be possible to obtain a 
better approximation of this value, and also to compute the coefficients of the 
asymptotic expansion. Since our intention was to compare the results obtained 
by the analysis of the solutions of the inner equation, we have not pursued in 
this direction. 

Now we compute U)m(z) ■ By definition (42) of Wj n (z) and (25) of the operator 

A, 



u> in (z) = ATre 2 ^ z Re(<d{z)-A{d z (j> s ){z)-d z (t> s {z)-^{Q){z))+0 



• g-2iriz z 2r + 2 e ~ 2 ~* iz 

(52) 

where we have used that by Theorem 3.4, (i(z) — d z <j> u (z) = 0(c 27Tlz z r+1 ). 

For symmetry reasons, we choose z — —ip with p E [2.25, 7]. We have used 
long double precision in C for calculating s ' u (z), d z 4> a ' u {z) . The strategy was 
suggested in [GL01] 

• First we compute the formal series 4>n up to order JV big enough. We know 
that the solutions ^> s ' u are close to </>jv if \z\ is big enough. Analogously 
for d z <jf< u . 

• We evaluate the formal series 4>m{z ± k) and <9 2 0at(z ± k) with k G N big 
enough. 

• Since both cff' u satisfy the inner equation, we obtain cff' u (z) and (j> s - u (z + l) 
recurrently. Analogously for d z (jf' u (z) and d z cff' u (z + 1) 

We have computed w; n (z) for the inner equations (41). Our results are given 
in the following picture, where we have added the scaling factor esc = 871 • 10~ 5 . 



). 



6.772 



6.7707 



v 



6.7695- 



6.7682- 
-7 



-log(u in (h)) for /i 
4og(aj in (/i) ■ esc) for f 2 



-6.05 



-5.1 



Imz 



-4.15 



-3.2 



-2.25 



We can observe that, on the one hand when Imz E [—3, —2] the theoretical 
error in (52) is big. On the other hand, when Imz G [—7,-6] the round- 
off errors (for /i) begin to be bigger than the theoretical error and hence the 
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computed values have noise. Nevertheless for values of Imz € [—5, —3] the 
computed values of u>i n (z) behave like a constant. More precisely, we have 
found u> in (z) = 871.6833 . . . for A 2 (</>) = -0 7 and w in (z) = 1.000832 . . . 10 5 for 
A 2 (</>) = — </> 3 which agree with the results for uj(h) given in (50) and (51). 

5 Formal solutions of the inner equation 

In this section we prove the existence of formal solutions of the inner equa- 
tion (27). The proof of the existence of formal solutions of equation (28) follows 
the same procedure. Hence, we skip it. 

We start by defining the spaces which these formal solutions belong to. For 
neN and r(n — 1) = 2, we define 

X r ={0(z)=]r^i|c fc ec}, (53) 

fc>i z 

the space of formal power series in x~ r without constant term, and, if n — 2m— 1, 
m > 2, that is, r — l/(m — 1), 

X** = = £ £ i.,^ <>.,■ }• (54) 

where [x] denotes the integer part of x, the space of formal power series in 
x~ r and logz, with the power of log z bounded by the power of x~ r , without 
constant term. 

We will say that 4> = Cfcr, with k g N, if and only if z kr 4> e C[[z~ r ]] is a 
power series with terms z~ 3 ' r for j > 0. We will also use Ok r ,j in X^ s , with 
k £ N, j G N U {0}, meaning that (f> = Okr.j implies that z kr (log z)~ : '4>(z) is a 
formal power series with terms of the form z~ k T log J z, with k' > and j' > —j 
such that j' < whenever k' — 0. We keep both notations in order to emphasize 
that Okr is a series without logarithms, while Ok r ,o is a series whose leading 
term does not have logarithms. 

We collect several properties of these spaces in the following lemma, whose 
proof is straightforward. 

Lemma 5.1. Let n > 2, r = 2/(n— 1) and g an analytic function around the 
origin with g(y) = Ay + 0(y for some I £ N U {0}. The spaces X r , for n 
even, and X], os , for n odd, have the following properties: 

1. X r and X l ° s are invariant by the formal differential operator -§^s ■ Fur- 
thermore, if <j) € Xl° s (resp. X r ), then ^4>(z) = z~ 2 ip(z) with tp g X l ° g 
(resp. X r ). 

2. If 4>(z) = az~ r + 4>(z), with 4> = 2r ,j, < j < [l/(m - 1)], resp. 2r , 
then g(az~ r + <fi(z)) = Aa e z~ tr + <p(z), with ip = Oti+xy j, resp. Ott+\\ r . 

Moreover, in the case n = 2m — 1, X)° g is also invariant by translation, that 
is, if 4>(z) £ Xl° s , then <f>(z — c) g X}° s , for any c g C. In the case n = 2m, if 
4> g X r , then <p(z — c) g X r / 2 . 
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We recall the function g(y) = y n — G(y). We remark that, since the operator 
A 2 can be written formally as 

A**M = 4sinh 2 m = (*L + 1^ + • • •) (55) 

(1) in Lemma 5.1 implies that the inner equation (27) is well defined in X r and 
X l r °s. We introduce 

£(</>) =A 2 (0) -.9(0). (56) 

It is clear that e(</>)(z) = z~ 2 e(z) with e e X| og (resp. X r ). 
Next lemma follows directly from the definition of X l ° E . 

Lemma 5.2. Let n = 2m — 1, r = l/(m — 1) and <\> £ Xl° s and e(<fi) as 
in (56). // z 2 t{(j)){z) has no terms of order N or smaller in z~ r (that is, no 
terms of the form z~~ kr log J z, with 1 < k < N), then e(cf>) = C(jv+i)r+2,L> where 
L = [N/(m- 1)]. 

Definition 5.3. Let n > 2, N £ N, r = 2/(n - 1) and <p £ X r or X^. We 
will call truncated series of order N of <f> to 4>N having the form: 

1. If n is even, 

N 

4> N {z) = 

fe=i 

2. If n = 2m — 1 is odd, 

m-1 ^ N I 

k=l ' k=rn 0<j<\£^\ 

Along the proof of Proposition 3.2, we will need to compute several times 
the formal series g(4> + ip) — <?(</>), with different cf> and ip. The following lemma, 
which follows from the properties in Lemma 5.1, summarizes the result. 

Lemma 5.4. Let n > 2, r = 2/(n - 1), N > 2, N £ N and </> £ X r or X l °%. 
We define ipN = 4>n ~ 4>n-i where 4>n and 4>n-x are the truncated series of 
order N and N — 1 respectively. We have that 

1. If n is even, 

g(<f> N (z)) - g{<j) N -i{z)) = -n—^tjj N {z) + C(jv+i)r+2- 

2. If n — 2m — 1 is odd, writing L = [N/(m — 1)] , 

g((j) N (z)) - g{<p N _ 1 {z)) = -n—^-ip N (z) + (N+1)r+2 . L - 

The following proposition implies the existence of formal solution of the inner 
equation (27) and henceforth Proposition 3.2. 
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Proposition 5.5. Letn >2,r~ 2/(n— 1) andco be such that Cq -1 = — r(r+l). 
T7ie inner equation (27) admits a formal solution <p with z r (4>(z) — c$z~ r ) € X r , 
i/n is even, and z r ((f>(z) — cqz^ t ) € Xl° s if n is odd. 

Let N > 2 and <f>N be the truncated series defined as in Definition 5.3. 
Writing the truncation error of order N as 

e N := e(4> N ) = A 2 (4> N ) - g(4> N ), 

where e was defined by (56), we have that 

1. If n > 2 is even €n = 0(jv+i)r+2- 

2. Ifn = 2m-1>2 is odd and L = [N/(m - 1)], then 

ft) ifl<N<m-l,e N = (N+1)r+2 , 
(ii) if m < N , e N = 0( N+1 ) r+2:L - 

Proof. We deal first with 1). We prove the claim by induction over TV. We start 
by assuming N = 1. Let (f>i(z) = c$z~ r '. By (2) in Lemma 5.1 and using (55), 
we have that 

- fin rii 

e 1 (z) = A 2 (<M(z) - g(M*)) = r(r + + ^ + («+D- 

The claim for N = 1 follows from the facts that r = 2/(n — 1), which implies 
C(n+i)r = 2r +2, and c£ _1 = -r(r + 1). 

Now we assume the claim for N — 1, that is, there exist coefficients Ck, 
1 < k < N — 2 such that 4>n-i satisfies 

£AT_i(z) = e(^JV-l)(^) = + °(N+l)r+2- 

We look for <Pn{z) = (j>N-i(z) + cn-iz~ Nt satisfying the claim. We have that 

e N (z) = eAr_i(z) + A 2 i^zwf) ~ 9 (j>N-i(z) + ^Zwf) + 9{4>n-i{z)) 
By 1) of Lemma 5.4, 

n— 1 

g (^> N -l(z) + C ^wf) ~ 9{4>N-l{z)) = ^Nr + °{N+l)r+2- (57) 

Hence, using again (55), 

, . _ AtV-1 . , CN-l rn 
eNi - Z > ~ z Nr+2 + AN z Nr+2 + U (N+l)r+2, 

where the coefficient Ajv is 

X N = 7Vr(iVr + 1) + nc^" 1 = * (n - ^ ) (N + n). (58) 

Clearly, the claim follows if Ajy is different from 0, which is true since n is even 
and positive. 

Now we assume n — 2m — 1, m > 2. The induction process used in the 
previous case can performed provided that Ajy 7^ 0. This is true for N =/= m. 
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Hence, the claim holds for 1 < N < m — 1. Let <f> m _i(z) = cq/z t + ■ ■ ■ + 
c m -2 1 ' z^ m ^ r be the corresponding function. It satisfies, 

A _i 

e m -i{z) = e(0 m _i)(z) = ™ r+2 + {m+1)r+2 - (59) 

Now we consider the case N = m. Since A TO = 0, this case cannot be dealt 
as before. We need to include logarithms in the formal series. 
Notice that, from (55), 

(60) 

We look for <f> m = (j> m -i+ip m satisfying the claim, with Tp m (z) = c m -i,iz mr logz+ 
Cm-ifiZ~ mr ■ Hence we have that 

e m = e m _i + A 2 (ip m ) - g(4> m -i + i> m ) + g(4> m -i). 
From (60), we have that 

A 2 (^ m )(z) = mr ^ m J 2 + ^ ip m {z) - (2mr + l)c m - hl ^-^ + O mr+4A , (61) 
while, from 2) in Lemma 5.4, 

g(<j) m -i{z) + ip m (z)) - g{4> m -\{z)) = -n^-ipmiz) + (m+1)r+2 , L , (62) 

with L = [N/(m - 1)] = [m/(m - 1)]. 

Hence, substituting (61) and (62) into the expression for e m above, we obtain 

e m (z) = e m -i(z) + -^il) m (z) - (2mr + l)c m _ M * +2 + C(m+i)r+2,L, 

where the coefficient Ajv was introduced in (58) and, in fact, satisfies A m = 
0. Since e m -i(z) = A m ^z mr+2 + 0( m+ i )r+2 (see (59)), taking c m _ M = 
Am-i/(2mr + 1), we have that e m = 0( m+1 ) r+2 ,L- Notice that the coefficient 
c m _i.o is free. Hence, the claim is proven for 1 < N < m. 

Now proceeding by induction the result is proven. □ 

6 A solution of the inner equation 

The goal of this section is to prove the existence of a solution of the inner 
equation satisfying the properties stated in Theorem 3.3. 
For any 7, p > 0, we recall the complex domains 

D* jp = {z£C: |Imz| > - 7 Rcz + p], D™ p = -Dl /p . 

defined in (34) (see Figure 1). We also introduce the norms 

IMI2*, P = sup \z» V (z)\ 



Inner equation for area preserving maps 



25 



and the Banach spaces 

X% tP = W ■ D u ^ -> C such that |M|^ >p < -f-oo}. 
We also define the functional space 

Ki^.p = W ■ C such that Tp(z) := (log z)- k <p(z) e 

and, if there is no danger of confusion, we will simply denote them 

V >" u > s jylog y«|!i II || || II u,s 7-j n u > s 

"■v — *-v,-y,pi ^v,k ~~ ^u,k,y,pi I ' II" — II ' 11^,7, pi ~ 1,P' 

From now on we will denote by C a generic positive constant independent 
of 7, p, v. We state (without proof) the following lemma which will be used 
without mention along this section. 

Lemma 6.1. Let < v\,V2- For any tpi £ X Ul and (f2 £ X U2 , then 

(fi ■ (P2 G X Vl+I/a and \\tpi • ^alU+i*, < II V2IU • ll</>i|U- 
Also there exists C > smc/i £/ia£ i/0 < v\ <vi and ip G X v . 2 , then 
V G ^ and l^li^ < Cp-^-'^II^IU. 

As in previous section, we will denote by O u and Vi u a generic function 
belonging to X u and X°% respectively. 

Theorem 3.3 is rephrased in terms of the Banach spaces X„'* . in the follow- 
ing proposition: 

Proposition 6.2. Given 7 > 0, there exists po > such that for any p > po, 

the inner equation (27) (polynomial case) and (28) (trigonometric case) 

A 2 (0)=g(<f>) (63) 
have exactly two solutions </> u ' s of the form 

where <f>o is the truncated series of order n defined in (35), (36) and (37), 
depending on the case we are dealing with, and ip u ' s G X r ^_ 2 . 

The properties of <po we are interested in follow from Proposition 5.5. 

Corollary 6.3. Let us consider the remainder of order n: 

eo = e(0o)=A 2 ((/>o)-.g(0o) 

where 4>q is the truncated series of order n defined in (35), (36) and (37). 
For any 7 > there exists po big enough such that 

1. Ifn is even, (f>o = CqZ~ 7 +02r , in the polynomial case, and 4>q = | log(— rz~ 
O r , in the trigonometric one. 

2. If n = 2m — 1 is odd, for the polynomial case 4>o = cqz~ t + 02r + Omr,i- 
Notice that, since m > 2, in particular we also have that (f>o = c^z~ r + 
C2r,i- In the trigonometric case, we have that (f>o = |log(— rz~ 2 ) + O r + 
C(m-i)r,i which also implies that <f>Q = | log(— rz~ 2 ) + Ti i. 
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3. For any value of n we have that cq € X nr +2 ■ 

The proof of Proposition 6.2 is performed in two steps. In Section 6.1 we 
introduce a linear equation which is close to the first order variational equation 
of (63) with respect to <pQ. Such linear equation can be easily inverted in the 
adequate Banach spaces. Finally, in Section 6.2 we look for ?/> u ' s as a solution 
of a suitable fixed point equation. 

From now on we will only deal with the — u— case, being the — s— case 
analogous. For that reason we will skip — u— from our notation. 



6.1 The linearized inner equation 

We introduce the function 

H(z) = (l + z- 1 Y-2+{l-z- 1 Y (64) 

for both cases, the polynomial and the trigonometric one with £ defined in (39). 
In this section we are going to study the following linear homogeneous second 
order difference equation: 

A 2 (0)(z) = H(z)<p(z). (65) 

We recall that the wronskian of two solutions, </>i,</>2 of a linear difference 
equation is defined as: 



W{4>u4>2){z) 



</>i(z) <j> 2 {z) 
0i(z + l) 0a(2 + l) 



In addiction, on the one hand, equation (65) has the obvious solution 772(2) = z e 
and, on the other hand, it is a well known fact that r\\ = b ■ 772 is a solution of 
(65) if and only if 

Ab(z) 



m{z) ■ r? 2 (z + 1) z e {z + iy 

One can also deduce that W(r]i,rj2) = 1. 

We will need a right inverse of the linear operator A defined in appropriate 
Banach spaces. For this reason we introduce the formal operator 

A- 1 (h)(z) = J2 h (z-k) (66) 

k>l 

We emphasize that we are dealing with the unstable case. 

Lemma 6.4. Let a > 0. For any 7 > there exists po > such that, for any 
p > pa, A -1 : X a+ i n>p — > <^ Q ,7,p is a right inverse of the operator A defined 
m (25) with \\Ar l \\ < C. 

The proof of this lemma is straightforward and can be found in [Gcl99] . 
The first variational around <fio of the inner equation (63) is given by 

A 2 (</)) = A 2 (0) = Dg(<f> o )0 (67) 
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and we notice that 

+ DG{<Pq) polynomial case 



V9[<Po) ~ { -(n - l^"" 1 ) + DG(e*>)e* 



trigonometric case. 



By using the identities Cq 1 = — r(r + 1) and nr = r + 2, the fact that 
H(z) = (£— l)£z~ 2 + 03 and Corollary 6.3, the following result follows: 

Lemma 6.5. For any 7 > there exists po > big enough such that 

1. The function H(z) satisfies H — Dg((f>o) — - A, with A € X r +2, if n 7^ 3, 
and A <G X r °^ 2 i> */ n = 3, 

2. The function 772(2) = Z is a solution of equation (65). Consequently, the 
function r\\ defined by 

V1{Z) = + !))<' 

is afeo an independent solution with W (1)1,7)2) = 1. -By Lemma 6.4, T)i € 
Xl-i 

We notice that property (1) of Lemma 6.5 implies that the linear equa- 
tion (65) is a good approximation of the first order variational with respect to 
0o (67). 

Finally, as we will see in the lemma below, Lemma 6.5 allows us to invert 
the linear operator C((f))(z) = A 2 (<fi)(z) — H(z)<fi(z). 

Lemma 6.6. For any 7 > 0, there exists po > such that for any p > po, the 
operator £(<fi) = A 2 (</>) — H ■ <j> has right inverse C^ 1 : X a +2,-y,p — > X a , 7tP i if 
a > I — 1 and it has the expression 

C-\h) = r)i ■ A -1 (772 • h) - in ■ A _1 (?7i • h). (68) 

Moreover, ||£ _1 (/i)|| Q ,, 7 , p < C\\h\\ a+ 2,-f.p being C an independent constant of 

Proof. We will skip 7, p from the notation. On the one hand 771,772 are inde- 
pendent solutions of the homogeneous linear equation £(4>) = and hence, by 
the variation of constants method, we obtain formula (68). On the other hand, 
if g € X a + 2 with a > £ — 1, then, 772 • g <G X a+ 2-e and 771 • g <E X a+ i + \ and 
by Lemma 6.4, 771 • A _1 (?72 • g) <G X a and 7/2 • A _1 (t7i • g) € X a . The bound 
l|£ _1 (5)lla — C||.°I|q+2 is obtained by a direct application of Lemma 6.4. □ 



6.2 The fixed point equation 

In this section we are going to prove Proposition 6.2 about the existence and 
properties of solutions of the inner equation (27) (polynomial case) and (28) 
(trigonometric case) 

A 2 (0) = -g(<P) 

of the form tp = (po + ip, with 4>q given by (35) (n even), (36) (n odd, polynomial 
case) or (37) (n odd, trigonometric case). 
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We introduce 

e = -A 2 (0 o )+ 5 (0o), ft(V0 = W D 2 g(<j ) o + Xyj)(l-X)d\ (69) 

Jo 

and we note that if <f> — 4>q + ip is a solution of the inner equation then, by (1) 
of Lemma 6.5 ip has to satisfy the second order difference equation given by 

A 2 (ip)- H ■ip = e +A-iP + TZ(ip) (70) 

As we proved in Lemma 6.6, the linear operator C has a right-inverse in some 
adequate Banach spaces. Using it, we will obtain a solution of the equation (70) 
by using the fixed point equation given by 

i/> = F{ip) := Cr l {e Q ) + C~\A ■ V>) + C' 1 ° K{ip). (71) 

Proposition 6.7. Let 7 > 0. There exists p\ > big enough such that, for any 
P > Pii the fixed point equation (71) has a unique solution ip £ X r +2,-(,p- 

Proof. We first note that there exists po > such that £ _1 (eo) £ X r +2,-y,p since, 
by Corollary 6.3, eo <G AV+4, 7iPo if po is large enough. Let go = 2||£~ 1 (eo)||r+2,7,p 
Along the proof of this proposition we will denote by if a generic constant de- 
pending only on <j>o, po and 7 and we will omit the dependence on 7 and p in 
the Banach spaces and norms. 

Let ipi,ip2 € B(qo) C X r+2 - We start by bounding the difference H-F^i) — 
J r (i/'2)||r+2. By Lemma 6.5 we have that, taking v r = 0, if n 7^ 3, and v r = r/2, 
if n = 3, we have that A £ X r ^-2-v T ,'y,p 1 provided that p\ is large enough. 
Henceforth, if ip £ X r+ 2, A ■ ip £ X2 r -\-&-v T - Applying Lemmas 6.1 and 6.6 we 
can easily check that 

\\C-\A ■ (^ - V 2 ))||r+2 < Cpi T+V r\\A\\ r+ 2- Vr ■ Ul - V' 2 ||r+2. (72) 

Now we deal with 1Z(ip{) — lZ(ip2)- We recall that 7Z was defined in (69). We 
notice that 

ft(^)-ft(^)=(V?-^a) / D 2 g(<f >0 + \yj 1 )(l~\)d\ 

Jo 

+ ^2/ [^ 2 .g(0o + AV'i)-^ 2 5 (0o + AV'2)](l-A)dA. (73) 
Jo 

We first claim that, if A £ [0, 1] and z £ D 7 Pl with p\ big enough, 

\D 2 g(Mz) + XM*))\ < K \ z f~ 4 < K\z\~ 2+r (74) 

where I was defined in (39). Indeed, we deal first with the polynomial case. In 
this case, by definition (29) of g, there exists a constant K such that \g(y)\ < 
K\y\ n . Moreover since g is an analytic function, Cauchy's theorem implies that, 
if y £ D(f>/2) 

\D 2 g(y )\ < K\y \- 2 sup \g(y)\ < K\y \ n - 2 . (75) 
\y-ya\<\yo\/2 

Also, since ipi £ B(qq) C X r+ 2, there exist constants < K\ < K2 and p\ big 
enough, Xi|z|" r < \<p (z) + Xipi(z)\ < K 2 \z\~ r < g/2 for any A £ [0,1] and 
z £ D-y, Pl ■ Then, using nr = r + 2 and estimate (75) 

\D 2 g(Mz) + AVi(*))| < K\M*) + A^iWI"" 2 < K\z\ 2r \z\- rn = K\z\~ 2+r 
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which proves bound (74) in the polynomial case. The trigonometric case is easier 
since \g(y)\ < K\e y ^ n ~ 1 ^\ and henceforth, a standard Cauchy estimate leads to 
bound (74). Hence, if tpi,ip 2 G B(g ), 

\D 2 g(Mz) + A(lM*))) ' - ^(*))l < K\z\-^{z) - Mz)\- 

Now we claim that, for A € [0, 1] and ipi,ip 2 G B(qo), 

\[&g(M*)+KM*)))-DMM*)+HM*))M^^^ 

Indeed, since g is an analytic function, D 3 g is bounded in D(p) and henceforth, 
for any yi, y 2 G D(p), |L» 2 .g(yi) - D 2 g(y 2 )\ < K\y\ - y 2 \ and the claim is proved 
provided that p\ is large enough to ensure that for any z G D 7 Pl , ipi (z), ip 2 {z) G 
B(g). 

Finally by using the previous computations and formula (73), one obtains 
that lZ(ipi) — Ti{ip 2 ) G X r +& U <% , ( 2r +4)+ r +2 = <%V+6 C A" r +4 and moreover 

\\n{i> x ) - n^ 2 )\\ r+i < cipi- 2 !^! - ^||r+2. (76) 

Then, by Lemma 6.6, £ _1 (7?.(V>i) — 1Z(ip 2 )) G A? r+ 2 and moreover 

Hr-^^^l) - ft(V> 2 ))||r+2 < <?|p|- 2 ||^l - ^Hr+2- 

Using this bound, (72) and definition (71) of the operator F, one has that, if p\ 
is large enough and p > p\ 

ll-F(Vl) - H^)\\r+2 < CpT^H^ - ^ 2 ||r+2 < ^1 - ^ 2 ||r+2 

and hence J 7 is contractive (we recall that i — v r > 0). Moreover, if ^ G B{qq), 

||^WIIr+2 < 11^(0)11^+3 + ||^(0) - J-W||r+2 < ||e || r +2 + \U\\r+2 < Qo 

which ends the proof of the Proposition. □ 

7 The difference U — <jf 

By Proposition 6.2 the existence of two solutions cf> u ' s = (j>o + tp u ' s of the inner 
equation is proved. Let us write 6 = (f> n — (j> s and we also introduce the function 

E=- [ (1 - X)D 2 g{<f> s + A(0 U - s )) d\- (<t> n - s ). (77) 
Jo 

We recall that both <fi u,s are solutions of the same nonlinear difference equa- 
tion: 

A 2 (0) + G(0H -.<? WO . (78) 

Consequently, the function satisfies the linear difference equation 

A 2 (e) = (-Dg(<j) s )+E)-e. (79) 

Although we do not have a good representation of the difference = <fi u — s , by 
means of Proposition 6.2 we already know it is well defined and some not optimal 
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bounds for which allow us to define a new linear equation from which is 
also a solution. In conclusion, we will use = (j) u — (f> s both as a known function 
(to define E(z)) and as an unknown solution of the above linear equation. 

The goal of this section is to prove that any analytic solution of equation (79) 
satisfying adequate boundary condition has to be exponentially small, that is 
of C(e _27rlz ). In fact, as claimed in Theorem 3.4, we will provide an exact 
formula for 0. 

7.1 Notation 

Given p, 7 > 0, let us recall the complex domain 

E 1}P = D\ p n D* jp n{zeC: Imz < 0}\{z G C : |Rez| < 1, |Imz| < p + 7} 

defined in (38) (see Figure 2). 

For i/.keR, we also introduce the norms 

|MU 7 ,„= Sup \z"<p(z)\, |Mfi.y,„= SU P \z V Qogz)- k <p(z)\, 

Z £ Ej *Y - p Z £ Ej—^ ^ p 

and the Banach spaces 

yu,i,p = {<P ■ Ey^p — > C such that ||9s||i/, 7 ,p < +00}, 
yi° 6 kn , P = W ■■ E y , p -> C such that |M|£f i7 , p < 
If there is no danger of confusion we will simply denote 

y» = y^,t,p, 11 • 11 1* = 11 • iu,7,pi y v % = y v ^^,p^ II ■ IL,I = II ■ lL,I,7,p- 

Lemma 7.1. Let < v\, V2- For any f G y vi and g G y„ 2 , then f ■ g G y vi + V2 
and 

11/ ■ 3lk+^ 2 < ll/lk ' IMU- 
j4/so, i/iere exists a constant C such that if < z/i < v<i and f G y v2 

fey vi and \\f\\ Vl <Cp-^-^\\f\\ U2 . 

As in previous sections, we will denote by O v and O v .k a generic function 
belonging to y vn , p and y„°f i7jP respectively. 

7.2 A right inverse of the operator A(<p)(z) = (p(z+l) — (j>(z) 

In this section we are going to construct a right inverse of the linear operator 
A: 

A(0)(*) = 0(« + l)-#z), (80) 

defined on functions belonging to 3^°f with v, k G R. We will follow the results 
introduced in [FS01] (which provide an explicit formula for A -1 ) and we also 
give useful properties of this operator when it acts on 3^°pi % ■ 

We first notice that, since -E 7 , P is an open set, for any z G E~ l ^ p there exists 
a{z) such that {w G C : \z — w\ < 2a(z)} C E ltP . In consequence, the complex 
path 7 2 = 7^ V 7^ (see Figure 3) 

l 1 z {t) = {-i{p + 1 ){l-t) + t{z-a(z)), (6 [0,1)}, 

7|(t) = {.z-cr(z) + it, £G(-oo,0]} (81) 
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is contained into the complex set E ltP . 

Given h an analytic function and z £ Ej tP , we introduce the linear operators 

A- 1 (h)(z)= [ - . kiu ] du and Al 1 (h)(z)=[ ^ r du. 

(82) 

Proposition 7.2. Let v, k £ R and 7 > 0. We define the linear operator 



A~ 



Al 1 ijf 1/ < 



A+ 1 if v > 0. 

There exists po > smc/i £/iai, /or any p > po, 

1. if v ^ 0, A -1 : yj+i 7 p - > 3^i°f 7 p * s a riffW inverse of the operator A. 

U. if v = 0, A -1 : 3^j°® — > 3^o°fc+i 7 p * s a r *ff^ inverse of the operator A. 

Moreover in both cases, there exists a positive constant C such that || A _1 || < C. 

Proof. Along this proof we will denote by if a generic constant depending only 
on 7 and v. We will skip 7, po and p from our notation of the Banach spaces 
and norms. 

We fix v £ R, 7 > fulfilling the hypotheses of Proposition 7.2 and po < p big 
enough. Let h £ y°?i k , and we introduce tp = A _1 (/i). Our first observation is 
that (p is an analytic function defined in E~ hP . Indeed, for any o~q > we define 
the set 

VL Uo = {u £ £ 7i p : u - do £ S 7i p}. 

We emphasize that E ltP = U <7 . >ofi(r . Moreover we note that, if z £ O CT0 we can 
take a{z) = er in the expression (82) of <f(z). Henceforth, in order to deduce 
that if is an analytic function in f2 CTo , we only have to study the convergence of 

-du 



e =F2fi(u-z) — 1 
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To this end, we observe that \ e T^ih 2 z -^) - 1| > | c ±t27r| P+7 +im 2 | _ ^ and that 

H^(t))\ < C(z)\t\-»- 1 log k (\t\) for some function C{z). Therefore, if v < 
and t £ (—00, 0], 

<2C(z)\t\-»- 1 log k {\t\)c t2 ^ +Inlz \ 



| e 27Ti( 7 J- 2 ) _ J] 



and we are done for the case v < 0. The case > can be done analogously. 

Now we arc going to check that A -1 is a right inverse of the operator A. 
We take into account that, if z, z + 1 € E ljP 

and therefore, since the only singularity of Th(u) / '( e =F2wi(u-z) — 1) i s u = z and 
it is a simple pole with residue /i(z)/27ri, we have that both ip± are solution of 
A(ip) = h defined in the complex domain E ltP . Here we have proceed exactly 
as in [FS01]. 

It only remains to prove that tp = A _1 (/i) e provided h € y°?i k - We 
restrict ourselves to the complex domain i? 7 P+7 C E ltP defined by: 

E y!p , = D^ p , n D s l p , n{zeC: Imz < 0}. 

We notice that, if the following bounds are proved, 

\z v Qogz)- k <p(z)\<K\\h\\% tk zeE^ p+1 if i/^0 (83) 

\{\ogz)- k -\{z)\ < K\\h\\ l ° s k zeE 7 , p+1 if u = (84) 

the same statement holds for z e -E 7 . p . Indeed, assume that bounds (83) and 
(84) arc satisfied and let z G i? 7iP \-E 7 . p . 7 . We have two cases, Rez < and 
Rez > 0. On the one hand, if Re < 0, it is clear that z + 1 E S 7iP+7 and that 
ip(z) = ip(z + 1) — h{z). On the other hand, if Rez > 0, z — 1 G -E 7 , p + 7 and 
consequently tp(z) = <p(z — 1) + h(z). In any case, \(p(z)\ < \(p(z ± 1)| + |/i(z)|. 
Here we have used that A(<p) = h. Therefore, if v ^ 0, using bound (83), we 
obtain 

|^(z)|< A-ll/ill^dziiniog^l^ + lzl — ^logz^) <K\\h\t%\z\- v \\ogz\ k 

and the result is proved for v 7^ 0. Analogously we check the result for v = 0. 

The proof of bounds (83) and (84) is easy but it requires tedious compu- 
tations which will be omit here. Nevertheless, we point out that for any fixed 
z E E ltP+7 we can take <j(z) = 1/2 in the definition (81) of j z . □ 

7.3 Two independent solutions of the linear equation (79) 
We recall that O = <p u — s satisfies equation (79): 

A 2 = (-Dg(<f> s ) + E)Q. (85) 

The following lemma states the properties of E we will need. Its proof is com- 
pletely analogous to the one of bound (74) provided U — eff G y r +2,-y.p- 
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Lemma 7.3. Let 7 and p satisfying the conclusions of Proposition 6.2 and E 
be the function defined in (77). We have that E G 3^ r +6-£ l7lP . 

As we did in Section 6.1, we split 

-Dg{4> s ) +E = H + M 

where H was denned in (65), M G 3V+2, if n 7^ 3, and M G 3^+2 u if ^ = 3. 
We rewrite equation (85) as 

A 2 (9) - H Q = M e. 

A solution of the homogeneous equation A 2 (ip) = H ■ ip is 772 (z) = z l . The 
function 

m ( z ) = /a- 1 ( 1 ) g y t -! 

is another solution satisfying W (771, 772) = 1. 

By using these decomposition as well as Proposition 7.2 for the operator A -1 , 
we can obtain solutions of the non homogeneous linear equation A 2 (ip) — H -ip = 
h. 

Lemma 7.4. For any 7 > there exists po > large enough such that for any 
P > P01 the operator C((p) = A 2 (1^2) — H • <p has right inverse defined in E ltP : 

£-\h) = rn ■ A" 1 { m ■ h) - r) 2 ■ A" 1 ^ ■ h). (86) 

There exists C > such that for any a G K and h G y a +2,-y,p, we have: 

1. If a 7^ £ — 1 and a 7^ — £ then £ _1 (/i) G iVa,7,p a^c? ||'C _1 (^)||q,7.p < 

C||/l||a+2, 7 ,p. 

«. // eitter a = £ - 1 or a = -£, ^(h) G ^°!, 7 , p and ||£- 1 (/0llL° g i l7 ,p < 

C||ft.|| Q +2,7,p- 

Next lemma provides a fundamental system of solutions of the linear equa- 
tion (85). 

Lemma 7.5. Let 7 > 0. There exists po large enough such that for any p > po, 
the equation (85) has two independent solutions, fji and f)2, satisfying 

m(z) = d z cjf(z) + r)i(z), fj[ G y r +3,7,P> 

UZ)= r dl (2£-l) + rt {z) > $ €y &™ 

with v = mm{r — £, 1 — k = 0ifn^=3,k = lifn = 3 and de is defined 
in (39), . 

Proof. First we look for 771. By construction, is a solution of the variational 
equation A 2 ^ = —Dg(<fr s )ip, therefore, the equation that fy\ has to satisfy is 

A 2 ((p)-H-tp = M ■ ip + E ■ d z <j) s . (87) 

We look for i)\ by means of the fixed point equation 

ip = £ _1 (J5 • d z (f> s ) + C~ l (M ■ if). 
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We are interested in solutions belonging to 3V+3- It is enough to check that the 
norm of the linear operator Q : y r +z — > defined by G(tp) = C^ 1 (M ■ ip) 

is less than one. This fact follows from Lemma 7.4 together with the fact that 
E G y r +6-t, M G y r +2 if n ^ 3 and M G 3^+2,1 if " = 3 - 0ne easil y tncn 
deduces that 

f,\ = (Id - Q)- 1 {C-\E ■ d z ^)) G 3W 3 (89) 

is a solution of equation (87). 

Now we deal with the second solution of the equation (85). We observe that, 
since 771 is a solution, then the function 7)2 = b ■ f\\ is also a solution of the linear 
equation (85) satisfying W ^2) = 1, if and only if b satisfies 

A(b)(z 



fji{z+l)fjx(z)' 

By Proposition 6.2, fji = d z (jf+f]\ = d z (j>n + d z ^ B +r)\ where d z ip s , f)\ G 3V+3- 
Moreover, using the definitions of dp,£ in (39) and Corollary 6.3, it is a direct 
computation to check that b has to satisfy the linear equation 

2^—2 

A (W*) = 75^ + s{z) 

with S £ 3^-2^+3 if n — 2, S G 3^2£ +3 ,i if re = 3 and S G 3^-2^+2+r if n > 3. 
We take 

„2«-l 

&o(*) 



*tfi(2€ - 1) 



and we note that r 2 d 2 A(bo)(z) — z 21 2 + 0-2f+3- Henceforth, the difference 
b\ = b — 60 satisfies an equation of the form 

A(6i) = S(z) (90) 

with S G 3^- 2 £+3 if n = 2, § G ^f f+3)1 if n = 3 and 5 £ 3^-2£+2+r if n > 
3. Applying Proposition (7.2) one has that equation (90) has a solution 61 
belonging to J 7 - 2^+2 if n = 2, 3 7 !°|^ +2 1 if ft = 3 and if tt, > 3 and the 

result follows. □ 



7.4 A final formula for 6 = cf) n — (ff 

Since 9 is a solution of the linear homogenous difference equation (85), the 
general theory allows us to write it as 

Q(z) = pi(z)rn{z) + p 2 (z)m{z) (91) 

with 771,772 two independent solutions of (85) and pi,P2 1-periodic, analytic 
functions in Ey tP . Moreover, if ^(771, 772) = 1, the functions p\ and P2 are 
determined by 

Pl (z) = W(@, V 2)(z), p 2 {z) = -W(Q 1 n 1 )(z). (92) 

Lemma 7.6. Let 7, p > and 771,772 two independent solutions of the linear 
difference equation (85) satisfying that 1^(771,7/2) = 1 and that 771 G y r +i md 

m G y-i. 
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Then there exist coefficients p\,p\ (depending on 771,2 J such that 

Q(z) = miz^Pi^+iniz^P^*- (93) 

fe<0 fc<0 

Proof. We first point out that, we already know that 8 = 4> u — (f) s = -0 U — ■0 s € 
3^r+2,7.p provided that, by Theorem 3.3, ip u ' s € I n addition, if h G y v ^,p, 

then A(h) G y v +\,2-i.2p- Indeed, standard arguments can be used to prove that, 
if h G y^,-y. P then d z h G y v +i,2-y,2p (sec for instance [Bal06]). Therefore, if 

z,z + 1 g E 2l ,2 P , 

\h(z + 1) - h(z)\ < \\d z h\Ui J q Vz+ \u+i < * II^IU+ij^p+r 

Using the above property, that 6 G y r +2 and formula (92) for p\,p2, one has 
that pi G 3^1 and P2 G 3V+4- In particular, pi,£>2 — ► as Imz — > —00 and since 
they are 1-periodic: 

Pl (z) = 5>fe 2ffifc *, p2(«) = J2P2 e ^ ikZ 

fe<0 fe<0 

ant the lemma is proved. □ 

We recall that the existence of independent solutions of the linear differ- 
ence equation (85) satisfying the hypotheses of Lemma 7.6 is guaranteed by 
Lemma 7.5. Henceforth Lemma 7.6 applied to 771,772 already gives an expres- 
sion of which is exponentially small. Among other things, we have proved 
that there exist 7, p > such that 

\e 2 * iz z- e (r(z) ~ <P(z))\ <K ze E^ p . (94) 

Nevertheless we have not proved Theorem 3.4 yet. We need to look for more 
suitable independent solutions of (85) to apply Lemma 7.6. 

Corollary 7.7. Let 7 > 0. There exists po > big enough such that for any 
P > PO; equation (85) has two fundamental solutions of the form 

Ci(z) = d z q> s (z) + ( 1 1 (z), 

^) = ^T^W + ^). 
with with di defined in (39) and (,\ satisfying 

sup l^-V^CiCOl < +oo- 



In addition, <G y r -e,i, P if n > 3, G yi-e :1 , P ifn = 2 and Cl € y^-i^^ if 
n = 3. 

Proof. As in the proof of Lemma 7.5 we write Ci = Ci ~ d z cjf. We note that £| 
satisfies equation (87): 

A 2 (Ci 1 )(z) - H(zKl(z) = M(z)Cl(z) + E(z)d z <j?{z). 
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We write = e 2mz ^(z) and we notice that £ has to satisfy the equation 

A 2 (C)(z) - H(z)C(z) = M(z)C(z) + c 2 ™E(z)d z cj) s (z). (95) 

We introduce ifa(z) = e 2nlz E(z)d z (jf(z). We first claim that ip € y^-i- 
Indeed, we note that d z <fi s G 3^-i and we recall that 

E(z) = - [\l-\)D 2 g(^ + \(r-<t> s ))d\{4> u -<t> s )- (96) 
Jo 

The claim follows from the facts that, by (94), \z~ e 2mz (4> u — <jf)\ is bounded 
and moreover \D 2 g(<fr s + \(<fr u — 4> s )) \ < K\z\ l ~ A if z <G E ltP (which can be proved 
as in (74)). 

It is clear that, a particular solution £ of (95) is given by a solution of 

( = £- 1 (vo) + G(() 

being g{0=L~\M- C). 

First we observe that, by Lemma 7.4, the independent term £ _1 (<y9 ) G ZVi— ^. 
Secondly we check that (Id — Q) is invertible in 3^1-^- Let ijj G 3^i-f. Since 
M G y$+2 for any n > 2, we have that M ■ ip G 3^3+ and consequently, by 
Lemma 7.4, £7(V0 G 3^i+r_f and moreover 

liaWOIk-/ < p- r/2 \\gm\i+^ < p7 /2 cu\\i-i < 

This implies that the norm of the linear operator Q : 3^1— e — > yi—t is less than 
one and therefore Id — Q is invertible. To this end, we can write £ as 

C = {id-gy 1 (c- 1 (^)) 

and we deduce that £ e yi-e and ||CI|i-£ < 2||£ -1 (< / 9 )||i-£ which implies the 
result for £i. 

The existence and properties of C2 follow from the ones for 7)2 in Lemma 7.5. 

□ 
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